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ABSTRACT 

Rest-frame 0.48-1.1 fim emission line strengths and molecular gas mass (H2) 
upper limits for three luminous infrared sources - the hyperluminous infrared 
galaxiesQ IRAS F09105+4108g {z = 0.4417), IRAS F15307+3252 {z = 0.926), 
and the optically-selected QSO PG 1634+706 {z = 1.338) - are presented. 
Diagnostic emission-line ratios ([O III] A5007/H/3, [S II] AA6716, 6731/Ha, [N 
II] A6583/Ha, and [S III] AA9069+9532/Ha) indicate a Seyfert 2-like spectrum 
for both infrared galaxies, consistent with previously published work. The upper 
limits of molecular gas mass for all three sources are M(H2) < 1 — 3 x 10^^h~^ Mq 
(go = 0.5, Hq = lOOh km s~^ Mpc~^), less than that of the most gas-rich infrared 
galaxies in the local Universe. All three sources have L\^/L'qq ~ 1300 — 2000, 
and thus are the extragalactic sources with the most extreme Li^/L'qq values 
measured to date. Given the relatively warm far-infrared colors for all three 
objects, much of their infrared luminosity may emanate from a relatively modest 
amount of warm dust (e.g. Md ~ 10^"^ Mq, Td = 200 - 100 K) near the AGN. 
For F09105+4108 and F15307+3252, the implied circumnuclear covering factor 
of this dust is ~ 90%, while for PG 1634+706 the covering factor is only ~ 35%. 



^HyLIGs: > lO^^L© where L„. = L{8 — lOOO/im) and assuming Hq = 75 km 

s~^ Mpc~^ and go = (e.g. Sanders & Mirabel 1996). Additional infrared galaxy 
definitions used throughout this paper - Luminous Infrared Galaxies (LIGs): galaxies with 

= IO^^^^^'^^Lq, Ultraluminous Infrared Galaxies (ULIGs): galaxies with L^,. > lO^^L©. 
Throughout this paper, we adopt Hq = lOO/ikm s~^ Mpc~^ and go = 0.5. 

^The source prefix 'F' will be used throughout this paper as a shorthand to refer to 
sources listed in the IRAS Faint Source Catalog (FSC: Moshir et al. 1992). 
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Subject headings: galaxies: ISM — infrared: galaxies — ISM: molecules — radio 
lines: galaxies — galaxies: active 
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1. Introduction 

Studies of ultraluminous infrared galaxies (ULIGs) detected by IRAS have focussed on 
revealing the nature of the host galaxies and the source powering the enormous infrared 
luminosities. Three forms of spectroscopic evidence for both stellar and nonstellar energy 
sources exists: (1) Studies of ULIGs comparing the fluxes of hydrogen recombination hues 
with optical and near-infrared forbidden lines have shown that the dominant ionization 
process appears to differ from galaxy to galaxy. Emission-line diagnostics have been used 
to classify these galaxies as H II region- like galaxies (i.e., emission lines induced by hot 
OB stars), Low Ionization Nuclear Emission-line Regions galaxies (LINERs), or Seyferts 
(e.g., Sanders et al. 1988b; Kim et al. 1995; Kim 1995; Veilleux et al. 1995). (2) Spectral 
polarimetry of ULIGs often reveal highly polarized continuum and broad emission-line 
spectrum, indicating the presence of a buried AGN (e.g., Kay & Miller 1989; Hines 
1991; Hines 1994; Hines & Wills 1993; Jannuzi et al. 1994). (3) Millimeter- wave (CO) 
spectroscopy of ULIGs has shown them to be rich in molecular gas, indicating the presence 
of an enormous amount of material for fueling an active galactic nucleus (AGN) and star 
formation (e.g.. Young et al. 1984; Sanders and Mirabel 1985; Sanders et al. 1987; Scoville 
et al. 1986; Sanders, Scoville, Soifer 1991; Sanders 1991). In addition, morphological studies 
of ULIGs have revealed that they commonly have tidal distortions and nearby companions 
indicative of the merger of two or more galaxies (e.g. Sanders et al. 1988a,b). Sanders et 
al. 1988a,b have suggested that ULIGs may be the precursor to optically-selected quasars. 

The recent discovery of the high-redshift HyLIG IRAS 10214+4724 {z = 2.286: 
Rowan-Robinson et al. 1991) and the subsequent detection of molecular gas in the source 
(Brown & Vanden Bout 1991, 1992; Solomon, Downes, & Radford 1992a), as well as 
evidence that it harbors a buried AGN (e.g. Elston et al. 1994; Jannuzi et al. 1994), 
has sparked interest in searching for and studying more luminous IRAS galaxies at early 
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epochs of the universe. To date, there are two other pubhshed findings of 2; > 0.3 galaxies 
discovered by IRAS - F09105+4108 {z = 0.4417: Kleinmann k Keel 1987; Kleinmann et al. 
1988) and F15307+3252 {z = 0.926: Cutri et al. 1994), and surveys are underway to search 
for additional candidates. In addition, a significant fraction of moderate-to-high redshift 
QSOs were detected by IRAS, making observations that test the infrared galaxy-QSO 
evolutionary hypothesis at higher redshift (and higher luminosities) feasible. To date, 
2 local QSOs (Mrkl014: Sanders, Scoville, & Soifer 1988a; IZwl: Barvainis, AUoin, & 
Antonucci 1989), one moderate-redshift quasar (3C48 Bit z — 0.3695: Scoville et al. 1993; 
Wink et al. 1997), and 2 higher redshift QSOs (the gravitationally lensed H1413-M17 at 
z = 2.556: Barvainis et al. 1994; BR 1202-07 at ^ = 4.695: Ohta et al. 1996; Omont et al. 
1996) have been unambiguously detected in CO. 

In this paper, we present CO observations and near-infrared spectroscopy of the high 
redshift HyLIGs F09105+4108, F15307+3252, and PC 1634+706 (IRAS 16347+7037: e.g. 
Sanders et al. 1989), the latter of which is the most infrared luminous, optically-selected 
QSO detected by IRAS in the PC sample. The purpose of these observations is to extend 
the determination of molecular gas properties and dominant ionization mechanisms to 
distant sources with extreme Lir, as well as to establish whether their properties are 
consistent with the trends in such properties observed for local ULIGs and QSOs. 

This paper is divided into 5 sections. The near-infrared and millimeter observing 
procedures are discussed in §2. The data reduction methods and fine- flux calculations are 
summarized in §3. Section 4 begins with a brief discussion of the optical and infrared 
morphologies of the three sources, then the dominant ionization mechanism in the two 
IRAS galaxies, the interpretation of the Lir and L'qq of the 3 sources, and their spectral 
energy distributions arc examined. Section 5 summarizes the results. Throughout this 
paper, we adopt Hq — lOOh km s~^ Mpc~^ and Qq = 0.5. 
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2. Observations 

Tables 1-4 summarize basic information and a journal of observations for the 3 sources 
discussed in this paper. Observations with each of the telescopes used are discussed 
separately below. A more in depth discussion of the observing technique and the difficulties 
inherent in millimeter- wave and submillimeter-wave observations of faint sources can be 
found in Evans et al. (1996). 

2.1. Imaging 

All of the imaging data were obtained at the UH 2.2m Telescope on Mauna Kea. 
The near-infrared imaging data were taken with the UH QUick Infrared Camera (QUIRC: 
Hodapp et al. 1996), which consists of a 1024x1024 pixel HgCdTe Astronomical Wide 
Area Infrared Imaging (HAWAII) array. The R-band image of F15307-F3252 and the 
I-band image of PG 1634+706 were taken with the UH Tektronix 2048x2048 CCD camera 
(Wainscoat 1996). Finally, the I-band image PG 1634-1-706 was obtained with the Orbit 
Semiconductor 2048x2048 CCD camera (Wainscoat 1996), reimaged to f/10. With the 
exception of the R-band image of F15307-I-3252 which consisted of one long exposure, all 
observations were done by taken a series of dithered exposures. 

2.2. Near-Infrcired Spectroscopy 

2.2.1. UH 2.2m Telescope 

All spectroscopic observations on the UH2.2m telescope were made with the K-band 
spectrograph (KSPEC: Hodapp et al. 1994). KSPEC is a cross-dispersed echelle configured 
to cover the wavelength range 1.1-2.5 iim in 3 orders {J, H, and K) on a, 256x256 
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NICMOS-3 HgCdTe detector array. The wavelength resolutions (at 2.2 //m) for each of the 
observations are listed in Table 3. Additionally, the wavelength range 0.7-1.0 /xm is also 
dispersed in several orders on the array. However, the order crowding at 0.7-1.0 /xm is such 
that only spectra of point-like sources can effectively be extracted, and this wavelength 
range is also compromised due to the fact that KSPEC is only in focus at A > 0.85 //m. 
In addition to its spectral capabilities, KSPEC also provides simultaneous imaging of 
approximately one square arcminute of sky around the slit on a second NICMOS-3 array. 
For the 1994 December observing period, the new UH tip-tilt system (Jim et al. 1997) was 
also implemented, notably improving the tracking, offsetting, and seeing of the observations 
and data. 

Each observing night began with a series of flat fleld exposures in the closed dome 
with incandescent lights turned on, then off. Once the telescope was guiding with the 
source in the sht, a 180-second exposure was taken, followed by a ISO-second exposure 10 
arcseconds off-source. Each spectrum exposure was accompanied by a 140-second image 
exposure. For the March observations, the pattern for observing was source-sky-sourcc-sky, 
etc., but the pattern was changed to source-sky-sky-source-source-sky-sky-source, etc. for 
the December observations to minimize spurious features created by the changing sky 
conditions. Observations of standard A (V) stars near each source were also made to flux 
calibrate the source spectrum and to remove telluric lines, and observations of an argon 
lamp were made to wavelength cahbrate the KSPEC data. 

2.2.2. United Kingdom Infrared Telescope (UKIRT) 

Observations of F15307-I-3252 were made with the upgraded Cooled Grating 
Spectrometer (CGS4) on the 3.8m UKIRT in shared-risk time. CGS4 is a 1-5 /xm 2D 
grating spectrometer with a 256x256 InSb array. Because the Hq;-|-[N II] A6583 and [S 
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II] AA6716,6731 (hereafter [S II] A6724) emission from the source is redshifted to J-band, the 
75 hne mm~^ grating in second order was used. Samphng was done by moving the detector 
over 2 pixels in 6 steps in the wavelength direction. For the observations, a 1 pixel-wide 
(1.2 arcsecond) sht was used. 

Spectral data was taken in 120-second exposures, with 20 seconds of each exposure 
made at each detector position during the 2-pixel sampling. The slit was then moved 15 
arcseconds in the spatial axis direction of the array for the next 120 second exposure. Thus, 
the source appears on the array as a positive and negative spectrum after the first frame 
is subtracted from the second. The source-sky observation sequence was the same as that 
used in the 1994 December KSPEC observations. Observations of the standard A star 
HD 162208 were made before and after the observations of F15307-I-3252. A krypton lamp 
was used as the wavelength calibration source. 

We encountered a problem with CGS4; the source flux of the positive spectrum was 
more than of the negative spectrum flux, implying that the source was only partially in the 
slit in the second position. Alternatively, the flux of the positive A star spectrum was less 
than the negative spectrum flux. The sht also appeared to be tilted relative to the array as 
evident by the position of the OH sky lines. Because of this, only half the data had enough 
flux to be usable, and no photometry was possible. 

2.3. Submillimeter Spectroscopy 

2.3.1. James Clerk Maxwell Telescope (JCMT) 

Observations with the 15m JCMT Telescope were made using the 1mm spectral hne 
receiver (A2), together with the Digital Autocorrelation Spectrometer (DAS) in wide-band 
mode (750 MHz bandwidth). Because of excess noise near the edge of the passband, the 
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usable bandwidth was only ~700 MHz, which corresponds to a total velocity coverage of 
~900 km s-^ at 230 GHz. 

All observations were obtained using a nutating subreflector with a chop rate of ~1.25 
Hz. Data were stored as six minute scans, and a chopper wheel calibration was performed 
after every other scan. Pointing was monitored every few hours by observations of standard 
continuum sources and was estimated to be accurate to within a few arcseconds. During 
the course of the observations, data were also taken with the bandpass centered at velocity 
offsets ± 40 km from the velocity corresponding to the redshift of the optical emission 
lines. Shifting the velocities in such a manner minimizes any ripples inherent in the baseline. 

2.3.2. Instituto de Radio Astronomia Milimetrica (IRAM) 30m Telescope 

The IRAM telescope has the ability to perform observations using the 3mm, 2mm, and 
1mm receivers simultaneously, making use of the 2 available 512 MHz filterbanks and two 
512 MHz autocorrelators. For F09105+4108, the 2mm and 1mm receivers were tuned to the 
redshifted CO (2 1) and CO (3 — > 2) transitions, respectively. During the observations, 
pointing was monitored by observing the planets and standard continuum sources and was 
typically accurate to within 3 arcseconds. 

2.3.3. National Radio Astronomy Observatory (NRAO) 12m Telescope 

Observations of the redshifted CO (2 — » 1) emission in PC 1634+706 were made with 
the NRAO 12m Telescope configured with a dual polarization SIS spectral-line receiver 
and two 256 x 2 MHz channel filterbanks. Observations were obtained using a nutating 
subrefiector with a chop rate of ~1.25 Hz. Six minute scans were taken, and a calibration 
was done every other scan. The pointing was checked every few hours by observing Jupiter 
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and was estimated to be accurate to within a few arcseconds. Because the duration of the 
observation was equivalent to one source transit, no velocity shifts were made. 

3. Results 

3.1. Imaging Data 

The I-band data reduction was done with IRAF. The data reduction consisted of 
flatfielding the individual images, scaling each image to its median value, and shifting and 
averaging the images. 

The K'-band data reduction was done similar to the I-band data reduction, except that 
after the median level of each individual frame was subtracted, the individual frames were 
averaged together (without spatial shifting) using a minmax averaging routine in order to 
create a "sky" image. This sky image was subtracted from the individual frames before 
they were shifted and averaged. 

Because the R-band image of F15307-I-3252 consisted of one integration, the image was 
simply flatfielded and scaled to its median value. Cosmic rays were then removed from the 
area around the galaxy. With the exception of the K'-band image of PG 1634-1-706, which 
shows only the unresolved point-like QSO, the final reduced versions of all the images are 
shown in Figures 1-3. 

3.2. Neeir-Infrared Spectra 

The KSPEC data reduction was done with IRAF. The sky frames were first subtracted 
from the source frames. The sky subtracted frames were then averaged and divided by the 
fiat fields. Because there were only a few bad pixels and none were on the spectral area 
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containing the source spectrum, the bad pixels were individually set to zero before the 
spectral orders were extracted with the APALL package. The extracted spectra were then 
wavelength calibrated and divided by the standard star spectrum, reduced in the same 
fashion as the source spectrum, to remove any instrumental effects and atmospheric lines 
(note that the stellar absorption hues in the standard star were removed using the IRAF 
task SPLOT before the division was done). Because of the narrow width of the slit, no 
photometry was determined for any of the KSPEC observations. Finally the spectrum was 
multiplied by a Planck blackbody spectrum with the same temperature as the standard 
star. 

Table 5 summarizes the emission-hne properties determined for F09105+4108 and 
F15307+3252. The emission-line fluxes, /(A), line widths, AfpwHM, and root-mean-square 
flux densities, 5'rms) were determined using the IRAF task SPLOT. The signal-to-noise 
ratios were determined by using the expression /(A)/Ai>FWHM'S'rms- 

The CGS4 data were reduced in a same manner as the KSPEC data, except that the 
flat flelding and masking of bad pixels was done automatically after every observation. 

Figure 4 shows the spectra of F09105-F4108 obtained with KSPEC. Six of the eight 
emission hues detected are at wavelengths corresponding to a source redshift oi z — 0.4417, 
identical to the observer-frame optical redshift determined by Kleinmann et al. (1988). The 
discrepancy in the [O III] AA4959,5007 doublet redshift {z = 0.4429) is due to wavelength 
calibration inaccuracies (the overlapping of the last few orders made it difficult to determine 
which argon line was in which order). 

Figure 5 shows the spectra of F15307-I-3252 obtained with KSPEC and the J-band 
spectrum obtained with CGS4. The KSPEC spectral fines imply a source redshift of 
z = 0.926, and the CGS4 spectra imply a redshift oi z = 0.927, both consistent with the 
redshift of 0.926 determined by Cutri et al. (1994). Note that both F09 105+4 108 and 
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F 15307+3252 show strong forbidden-line emission. Further, the Hq;+[N II] AA6548 + 6583 
and [S II] A6724 hnes widths of F09105+4108 (Table 5) are consistent with those measured 
from median resolution spectroscopy by Liu, Graham, & Wright (1996). The interpretation 
of the line emission will be discussed in §4.2. 

To determine the fluxes of the Ha and [N II] AA6548+6583 hnes of F09105+4108 and 
F15307+3252, the line profiles were fit using the IRAF package SPECFIT. Initial guesses 
to the fiuxes and line widths were entered into the program, and the line fiux ratios and 
wavelength ratio of [N II] A6583 to [N II] A6548 were fixed to 3 and 1.005345 (i.e., the ratio 
6583/6548), respectively. In addition, the line widths of the two [N II] lines were fixed to 
be equal. Iterations were made until a convincing fit was achieved. The best fit models are 
fisted in Table 5. 

Figure 6 is a short exposure of the redshifted [O III] AA4959,5007 and Ha fines of 
PG 1634+706. The redshifts determined from the Avfwhm = 4600 km s~^ Ha line and the 
[O HI] A5007 lines are z = 1.339 and 1.337, respectively. The systemic redshift determined 
by Tytler & Fan (1992) is 1.3371. 

3.3. CO Spectra 

Table 6 summarizes the CO line data for the 3 sources. In the case of the IRAM 
30m observations, where data was taken at several widely spaced velocity offsets, the data 
are presented in separate spectral blocks. Separating the data in this manner avoids the 
creation of direct current (dc) offset-induced ghost features (see Evans et al. 1996) and 
provides data blocks with relatively constant signal-to-noise across each spectrum. 

Figures 7-9 show the spectra of F09105+4108, F15307+3252, and PG 1634+706, 
respectively. For each source, individual scans were examined to check for sinusoids in 
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the baselines and other irregularities. After bad scans were removed, the rest of the scans 
were averaged and linear baselines were subtracted. For the JCMT 15m data, where the 
sources were observed with several small velocity offsets, only the overlapping regions of 
the scans were averaged. Finally, the spectra were smoothed to 30-50 km s~^. No emission 
lines were detected. More sensitive upper limits of the CO emission in PG 1634+706 and 
F15307+3252 have also been recently reported by Barvainis et al. (1998) and Yun, Scoville, 
& Evans (1998), respectively (see Table 6). 



3.3.1. CO Luminosity and H2 Mass Limits 

Table 6 lists the root-mean-square temperatures, Trms, achieved with a velocity 
resolution, Av^-es, for the 3 sources observed. To calculate the upper limit of the CO line 
intensity, Jco, we assume that the CO lines have a FWHM velocity, AfpwHM ~ 250 km s~^, 
similar to CO lines seen in local ULIGs (e.g. Sanders et al. 1991) and in IRAS 10214+4724 
(e.g. Solomon et al. 1992a), Using the observed Trms and At>res from Table 6, the 3a limit 
on Ico is given by 

Ico < ^= [K km s J. (1) 

res 



Multiplying Iqo by the Kelvin-to-Jansky conversion for a point source^ gives the 
corresponding limit for SqoAv [Jy km s~^]. 

The luminosity distance for a source at a given redshift, z, is. 



= cH,\,^ zqo + (go - 1) j2qoZ + 1-1 [h-' Mpc], (2) 



^ Kelvin-to-Jansky conversion factors: JCMT 15m- 15.6 Jy K^^ (Matthews 1992), IRAM 
30m - 4.6 to 4.9 Jy K'^ (Guelin, Kramer, & Wild 1995), NRAO 12m - 25.2 Jy K^^ (R 
Jewell, private communication) 
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where Qq is the deceleration parameter (= 0.5 for a critical density universe) and Hq — lOO/i 
(km s~^ Mpc"^). Given the estimated upper hmit of we can calculate the CO 

luminosity of a source at redshift z, 

^ (2^) *^coA^D^(l + z)-\ (3) 

where c is the speed of light, k is the Boltzmann constant, and z/obs is the observed frequency. 
In terms of useful units, L'qq for the C0( J J — 1) rotational transition can be written as. 

Table 6 lists the L'qq upper limits of the 3 sources observed. We have assumed that 
CO rotational transitions up to CO (4 3) have approximately the same brightness 
temperature. This assumption is based on previously published multi-transition CO 
analyses of the HyLIG IRAS 10214+4724 (e.g. Solomon et al. 1992a). For PC 1634+706, 
we have assumed that the CO (5 4) transition has approximately the same brightness 
temperature as the lower transitions, as observed in the Cloverleaf Quasar (Barvainis 1996). 

The upper limits on H2 masses implied by the observed limits on L'qq can be 
determined by making the reasonable assumption that the CO emission is optically thick 
and thermalized, and that it originates in gravitationally bound molecular clouds. For 
molecular gas in gravitationally bound clouds, the ratio of the H2 mass and the CO 



luminosity is given by a = M(H2)/-Lco \/n{'H2)/Tb Mq (K km s~^ pc^)~^, where n(H2) 
and Tb are the density and brightness temperature for the appropriate CO transition (e.g. 
Sanders & Scoville 1987; Solomon et al. 1992a). Multi-transition CO surveys of molecular 
clouds in the Milky Way (e.g. Sanders et al. 1993), and in nearby starburst galaxies (e.g. 
Glisten et al. 1993) have shown that hotter clouds tend to be denser such that the density 
and temperature dependencies tend to cancel each other. The variation in the value of a 
is less than a factor of 2 for a wide range of gas kinetic temperature, gas densities, and 
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CO abundances. We adopt a value of 4 Mq (K km pc^)~^ for a, which is similar 
to the value determined for the bulk of the molecular gas in the disk of the Milky Way 
(c.f. Scoville & Sanders 1987), as well as the value determined from a multi-transition CO 
analysis of IRAS 10214+4724 (Solomon et al. 1992a). Table 6 lists the upper limits on 
the the molecular gas mass of F09105+4108, F15307+3252, and PG 1634+706 determined 
using this value of a. 

4. Discussion 

As stated in the introduction, spectroscopy has been a powerful tool for determining 
the dominant ionization source in infrared-luminous galaxies. The environment and 
morphology of the host galaxy are an essential supplement to the spectroscopy, often 
providing information on global spatial scales that have direct bearing on the nuclear 
activity. We briefly discuss the imaging data, then dedicate the bulk of the remaining 
discussion to the spectroscopic data. 

4.1. Morphologies and Environment 

Figures 1-3 show images of the three sources taken at wavelengths longward of 
V-band. Most local {z < 0.3) ultraluminous infrared galaxies arc observed to have distorted 
morphologies, multiple nuclei, and close companions indicative of the merger of two or 
more galaxies (e.g. Sanders & Mirabel 1996 and references therein; Surace et al. 1997). 
F09105+4108 clearly shows evidence for companions (Figure 1). Indeed, unlike most 
luminous infrared galaxies (LIGs), F09105+4108 is clearly a cD galaxy. Both the I and 
K'-band images of the galaxy show evidence of "clumps" within the envelope of the cD 
galaxy, which may be galaxies stripped by tidal or dynamical friction. Such stripping may 
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be the source of dust for F09105+4108. In addition, F09105+4108 is known to have radio 
jets (Hines & Wills 1993), indicating the presence of a buried AGN. 

While the R-band image of F15307+3252 (Figure 2a) shows no direct evidence of 
distorted morphologies or close companions, the K'-band image Figure 2b) shows a faint 
companion to the southwest of the main galaxy. Soifer et al. (1996) and Liu, Graham, & 
Wright (1996) also present near-infrared images showing this companion, with the latter 
group showing evidence for a second companion. No evidence is seen for tidal features in 
Figure 2, or in any other images published to date. 

The QSO PG 1634+706 (Figure 3), the most distant of the three sources, shows no 
evidence of distortions or close companions. There appear to be approximately 30 faint 
objects in the vicinity of the QSO with mi — 23.5 — 25.5, which is consistent with the 
density of faint field objects expected within the given magnitude range (Pozzetti et al. 
1998). Thus, there is also no apparent evidence with the present data set that the QSO is 
imbedded in a cluster. 

4.2. Emission-Line Diagnostics 

Figure 10 shows 4 emission-line diagnostic diagrams commonly used to distinguish 
galaxies with Seyfert/LINER-hke properties from galaxies with H II region-hke properties. 
The usefulness of these diagrams can be understood by the following emission-hne 
characteristics: Emission from [S II] A6724 and [N II] AA6548-I-6583 can emanate from 
ionized hydrogen regions, and [S II] A6724 and [O I] A6300 emission can emanate from 
semi-ionized regions where coUisional ionization is significant. Enhancement of these lines 
occurs in AGNs because, unlike H II regions, they have extended partially ionized zones 
created by an excess of X-rays (the absorption cross sections of neutral hydrogen, helium 
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and all ions are small for X-rays, thus X-rays tend to escape the ionized region before 
interacting: Veilleux & Osterbrock 1987). [O III] A5007 is a high ionization line photoionized 
by UV photons, and thus tends to be strong in Seyfert galaxies. [S III] AA9069-I-9532 is 
not as sensitive to the level of ionization as [O III] A5007, most hkely due to its lower 
ionization potential (Osterbrock, Tran, & Veilleux 1992). However, [S III] AA9069+9532 and 
[S II] A6724 together do provide some degree of separation between the 3 classes of active 
galaxies (Figure lOd). 

Emission-line ratios of F09105+4108 and F15307+3252 are plotted in the diagnostic 
diagrams. For F09105+4108, the [O III] A5007/ H/? and [O I] A6300 / Ha ratios have 
been determined from Kleinmann et al. (1988). Because their [O HI] A5007/ H/3 ratio is 
being used in combination with our KSPEC measurements, it is worth cautioning that [O 
III] A5007/ H/3 may vary as a function of the aperture size used to observe F09105-I-4108. 
In all 4 diagrams, both galaxies lie in the Seyfert region, indicating that the dominant 
source of ionization is a hard, nonthermal continuum source. The analysis of these data 
are consistent with the detection of broad emission lines in polarized light in all three 
sources (Hines & Wills 1993; Hines et al. 1995; Jannuzi et al. 1994) and previous optical 
emission-line diagnostics of both galaxies done with low resolution (Soifer et al. 1994, 1996) 
and median resolution (Kleinmann et al. 1988; Liu, Graham, & Wright 1996) spectroscopy. 

Table 7 is a comparison of the observed [O HI] A5007 / Ha, [S II] A6724 / Ha, and [N 
II] A6548+6583 / Ha ratios of F09105+4108 and F15307+3252 with the sample of Bright 
IRAS galaxies and warm IRAS galaxies (which have been divided into Seyfert 2, LINER, 
and H II region classifications). For F09105+4108, the ratio [O HI] A5007 / Ha has been 
determined from Kleinmann et al. (1988), who calculate [O HI] A5007 / H/3 ~ 12 and 
a negligible extinction (note that the [O III] A5007 / Ha ratio calculated from our data 
indicate that half of the [O HI] A5007 flux is lost via defocussing; see §2.2.1). Despite the 
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obvious scatter in the [O III] A5007 / Ha ratio due to dust, there is a notable separation 
between the ratio for Seyferts and for LINERs and H II region galaxies (i.e., the separation 
is due mainly to the ionization parameter). While F09105+4108 has a [O III] A5007 / Ho; 
ratio 2 standard deviations above the average ratio for low-redshift Seyfert 2 galaxies, 
F15307+3252 have observed [O III] A5007 / Ha ratio higher than the average Seyfert 2 
ratio, but still within the standard deviation of the Seyfert 2 ratios. Inversely, the [S 
II] A6724 / Ha and [N II] A6548+6583 / Ha ratios are lower than the average values of the 
Seyferts, but still within the standard deviation of the Seyfert 2 ratios. This may imply 
that the hyperluminous galaxies have smaller semi-ionized regions or lower metallicities 
than the average, local IRAS Seyferts. 

4.3. Lir versus L'qq 

Given the evidence presented here and elsewhere for AGNs in both F09105+4108 
and F15307+3252, we consider the infrared and CO luminosities of these 2 sources and 
PG 1634+706 relative to LIGs, ULIGs and QSOs detected by IRAS, but having lower 
infrared luminosities. Figure 11a is a plot of Lqq vs. Li^.; the upper limits on the L'qq of 
F09105+4108, F15307+3252, and PG 1634+706 are less than the L'^q of the most gas-rich 
galaxies observed locally, perhaps indicating that the upward trend of data points in the 
diagram levels off at higher Lir, or disperses. Figure lib is a plot of L[t-/L'qq vs. Li^; the 
lower limits on the L^^/Lqq of the 3 sources are consistent with the increasing values of 
Lir/L'co ^ ^ function of increasing Lir of the sample of LIGs and ULIGs. One possible 
explanation for the data trends observed in Figure 11 is that Lir becomes increasingly 
dominated by reprocessed nonstellar light at high Lir. From the spectral energy distributions 
plotted in Figure 12, it is clear that all 3 sources have i/L,^(60//m)/i/L,^(100//m) ratios higher 
than nearby, less luminous ULIGs that have been detected in CO such as Arp 220 and 
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Mrk 1014 (e.g. Sanders et al. 1988b). The mass of "warm" dust required to produce the 
observed infrared luminosity can be computed by fitting the rest-frame far-infrared (40-100 
yum) emission with a single-temperature dust model using dust emissivities from Draine & 
Lee (1984). The dust temperatures implied by the far-infrared SEDs in Figure 12 are in the 
range 100-200 K, and the corresponding dust masses, = (Lir/10^ L0)(4OK/Td)^, are 
in the range 1O'^-1O^M0. Assuming, as is found for ULlGs like Arp 220 and Mrk 1014, 
that the bulk of the dust is heated to the computed "warm" dust temperature, and that 
the gas/dust ratio is similar to what is found for other ULIGs such as Arp 220 (i.e. ~200: 
Sanders, Scoville, & Soifer 1991), then the imphed range of H2 gas mass is 10^-^ — 10^-^ Mq. 
Thus it is possible that the observed infrared luminosities are produced by a relatively 
modest mass of gas and dust (e.g. comparable to or less than the total molecular gas 
mass of the Milky Way) heated to temperatures 100-200 K, perhaps by the central AGN. 
Such small amounts of molecular gas would also explain our non detections of all three 
sources in CO (see Table 6). Further, while it is clear from the observed narrow emission 
lines and large infrared excess (Lir/Lboi > 0.9) of F09105+4108 and F15307+3252 that 
the circumnuclear covering factor of the dust is almost unity, the broad emission hues and 
smaller infrared excess (Lir/Lboi ~ 0.35) of PG 1634-1-706 indicate that the QSO has a low 
covering factor. 

It is curious that CO surveys of galaxies at cosmological distances have yielded no 
evidence of galaxies with L'qq much larger than that of the most CO luminous galaxies at 
lower redshift {z ^ 0.2) plotted in Figure 11 (i.e., Lqq corresponding to ~ 3-4x10^° Mq 
of molecular gas). To date, there exists little evidence that the H2 content of active galaxies 
evolves with redshift (^ 4). The extremely high H2 masses originally computed for 
IRAS 10214-1-4724 and H 1413-1-117 were in large part due to gravitational lensing of these 
sources by foreground objects (IRAS 10214+4724: Elston et al. 1994; Soifer et al. 1995; 
Trentham 1995; Graham & Liu 1995; Broadhurst & Lehar 1995; Serjeant et al. 1995; 
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Close et al. 1995; Eisenhardt et al. 1995; Downes, Solomon, & Radford 1995; H 1413+117: 
Barvainis et al. 1994), and CO surveys of 1 < 2; < 4 high redshift powerful radio galaxies 
(HzPRGs: Evans et al. 1996; Downes et al. 1996) and z ~ 4 luminous high redshift QSOs 
(Barvainis & Antonucci 1996) were unsuccessful at detecting sources with more molecular 
gas than the most gas-rich galaxies observed locally (i.e. based on 3o" upper M(H2) limits 
of the HzPRGs and high redshift QSOs). The reported detection of CO emission in the 
z = 4.7 QSO BR 1202-07 (M(H2) ~ 6 x IO^^/i-^Mq) may be the first indication that 
M{B.2) increases beyond z ~ 4. However, reports of shear due to gravitational lensing in 
the vicinity of BR 1202-07 may mean that the CO luminosity is being amplified (e.g. see 
discussion in Omont et al. 1996b). Further, there have been recent, unconfirmed detections 
of CO emission in the z = 2.4 radio galaxy 53W002 (Scoville et al. 1997), the z = 4.5 QSO 
BR 1335-07 (Guilloteau et al. 1998), and the gravitationally lensed, z = 2.6 quasar MG 
0414+0534 (Barvainis et al. 1998), but all have molecular gas masses comparable to the 
most gas-rich, local galaxies. [] The inferred molecular gas masses and upper limits of high 
redshift galaxies surveyed to date may be evidence that either galaxies with more than ~ 
3-4x10^'^ h^"^ Mq of molecular gas rapidly turn the gas into stars, and/or that galaxies form 
in a hierarchical fashion. It may also be the case that only high-redshift galaxies with a 
significant nonstellar contribution to their infrared luminosity will be observed to have high 
Lij.. A better understanding of the relationship between extreme L[j. and H2 mass awaits the 
detection and observation of a larger sample of HyLIGs, especially those having relatively 
cool IRAS colors, i.e. relatively low values of z/L,^(60yum)/z/Lj,(100/im). 



''This statement assumes that the observed CO emission in MG 0414+0534 is 
gravitationally lensed by at least three times its intrinsic luminosity. 
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5. Summciry 

We have presented new millimeterwave C0(J+1,J) observations and rest-frame 
0.48-1.1 spectroscopy of the HyLIGs IRAS F09105+4108, IRAS F15307+3252, and 
PG 1634+706. The following conclusions are drawn: 

(1) The emission line ratios [O III] A5007/H/3, [O I] A6300/Ha, [N II] A6583/Ha, [S 
II] A6724/Ha, [S III] AA9069+9532/Ha of F09105+4108 and F15307+3252 indicate that 
both have Seyfert 2 emission-line spectra, consistent with previously published work. 

{2) Our upper limits on the molecular gas mass of the 3 HyLIGs are comparable to or 
less than the H2 mass of the most gas-rich galaxies observed locally. 

{3) The relatively "warm" far-infrared (60/xm/lOO/xm) colors (implying dust 
temperatures — 100 — 200 K) for all three sources are consistent with the idea that 
the bulk of their infrared emission emanates from a relatively modest amount of warm 
circumnuclear dust, i.e. = 10^ — 10^ Mq, corresponding to a molecular gas mass, 
M(H2) ~ 10^-^ — 10^-^ Mq. This would be consistent with our measured upper limits for 
L'qq and the corresponding extreme values of Li^/L'^Q. 

(4) While the narrow emission hues and Lir/Lboi > 0.90 of F09105-F4108 and 
F15307-I-3252 indicate a circumnuclear dust covering factor of almost unity, the broad 
emission lines and L„/Li,o\ ~ 0.35 of PG 1634-1-706 indicate that most of the AGN is 
exposed. 

We thank the staffs of the JCMT 15m telescope, the IRAM 30m telescope, the NRAO 
12m telescope, the UH 2.2m telescope, and the United Kingdom Infrared Telescope for their 
generous support during our observations. A.S.E. thanks N. Trentham, E. Egami, J. Hora, 
D. Kim, T. Greene, P. Hall, and J. Goldader for useful discussions and assistance during the 
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part of the K'-band and the R-band image data, respectively, of F15307+3252, and to J. 
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We also thank K. Teramura for the final preparation of figures show in this paper, and 
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NASA grant NAG5-3042. D.B.S. was supported in part by NASA grant NAG5-3370. This 
research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory. 



-25- 



Table 1. Source Parameters 



Name R.A. Dec. f„{X) (Jy) z Refs Li^'' 

IRAS B1950.0 12/im 25/im 60/im 100/im (Lq) 



F09105+4108'' 09:10:32.84 41:08:53.6 0.13±.03 0.33±.03 0.53±.04 <0.32 0.4417 1,3,4 7.8 x lO^^fe-a 

F15307+3252 15:30:44.63 32:52:50.6 <0.07 0.08±.07 0.26±.08 0.46±0.19 0.926 1,3,4 2.5 x lO^^h''^ 

F16348+7037 (PG 1634+706) 16:34:51.0 70:37:37.0 0.06±.01 0.13±.01 0.25±.03 <0.45 1.337 2,3,5 5.7 x 10^3^-2 

^Assuming qo = 0.5, Ho = 100/ikm 8"^ Mpc~l. 

''IRAS Point Source Catalog name in Kleinmann et al. (1988) was P09104+4109 

References. — Positions, fluxes, and redshifts: (1) Becker et al. (1995); (2) This paper (IRAS FSC); (3) This paper (IRAS SCANPI); (4) This 
paper (near-infrared emission lines); (5) Titler &c Fan (1992) 
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Table 2. Journal of Imaging Observations 



Source Telescope Detector Pixel Scale*^ Band Date Time 

arcsec / pixel ™/y min . 



IRAS F09105+4108 


UH 2.2m 


Tek 


0.14 


I 


03/97 


26 






QUIRC 


0.19 


K' 


04/96 


30 


IRAS F15307+3252 


UH 2.2m 


Tek 


0.22 


R 


04/94 


20 






QUIRC 


0.19 


K' 


04/96 


60 


PG 1634+706 


UH 2.2m 


Orbit 


0.14 


I 


06/97 


55 






QUIRC 


0.06 


K' 


05/97 


16 



^Imaging at the UH 2.2m Telescope is done with a focal ratio of f/10 or f/31. For 
the CCD observations, images taken at f/31 are binned 2x2 pixels. 
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Table 3. Journal of Near-Infrared Spectroscopic Observations 



Name 


Telescope 


Detector 


Slit Width 
II 


A 

AA 


Date 
mm/yy 


Time 
min 


IRAS F09105+4108 


UH 2.2m 


KSPEC 


1.0 


620 


12/94 


18 


IRAS F15307+3252 


UH 2.2m 


KSPEC 


0.8 


760 


03/94 


33 




UKIRT 


CGS4 


1.2 


1000 


05/95 


56 


PG 1634+706 


UH 2.2m 


KSPEC 


1.0 


760 


03/94 


3 
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Table 4. Journal of CO Observations 



Name line Telescope Date Time Tgy 

m/y hr K 



IRAS F09105+4108 C0(3 ^ 2) 

C0(3 ^ 2) 
C0(2 ^ 1) 

IRAS F15307+3252 C0(4 ^ 3) 

CO (4 ^ 3) 
C0(4 ^ 3) 

PG 1634+706 C0(5 ^ 4) 

C0(2 ^ 1) 



JCMT 15m 


04/94 


7.3 


429 


IRAM 30m 


06/94 


6.3 


724 


IRAM 30m 


06/94 


6.3 


530 


JCMT 15m 


03/94 


8.4 


380 


JCMT 15m 


04/94 


1.7 


400 


IRAM 30m 


06/94 


10.2 


590 


JCMT 15m 


04/94 


11.0 


452 


NRAO 12m 


05/94 


5.6 


160 



-29- 



Table 5. Near-Infrared Spectroscopy; Emission Line Properties 



IRAS F09105+4108 IRAS F15307+3252 



Line 


Instr. 


•^obs 


/(A)//([SII]) 


FWHM'' 


S/N 


-^obs 


/(A)//([SII]) 


FWHM» 


S/N 






(A) 




(km s-i) 




(A) 




(km s-i) 




H/3 


KSPEC 












< 1.4 


(1300) 




[O liI]A.1959 


KSPEC 


7153.1 


2.6±0.8 


740 


3.2 


9556.5 


5.5±0.4 


1540 


12 


\r\ TTTl \ r^nnT 
[\J iliJAoUUT 




7224.6 


9.0itl.O 


980 


8.9 


yD44.0 


16it0.4 


1280 


35 


[O iJAD3UO 






< 0.44 


(1000) 






< 1.0 


(1300) 




Ha + [N II] 


KSPEC 


9463.5*' 


6.8±0.2 


1710 


32 


12641.8 


9.8±0.5 


1990 


18 






9490.6<= 
















[N II]A6548<1 




9441.7±2.9 


0.62±0.09 


1160±150 




12618±5 


1.4±0.2 


1310±150 




HaA6563'' 




9460.9±1.9 


4.4±0.4 


1200±100 




12636±3 


5.2±0.5 


1300±150 




[N II]A6583d 




9492.2±2.9 


1.8±0.3 


1160±150 




12686±5 


4.1±0.4 


1310±150 




[S II]A6724 


KSPEC 


9693.9 


1.0±0.2 


1120 


6.2 


12960.3 


1.0±0.2 


850 


4.4 


[S III]A9069 


KSPEC 


13076.9 


0.74±0.08 


900 


8.7 


17467.1 


1.1±0.3 


630 


3.5 


[S III]A9532 


KSPEC 


13174.3 


2.0±0.08 


930 


20 










He IA10830 


KSPEC 


15615.5 


2.1±0.08 


1070 


27 


20873.8 


4.8±0.3 


1540 


16 


[O I]A6300 


CGS4 










12142.8 


< 1.6 


(1300) 




Ho + [N II] 


CGS4 










12650.4 


7.1±0.7 


1860 


10 


[S II]A6724 


CGS4 










12959.6 


1.0±0.3 


1130 


3.2 



^Uncorrected for instrumental broadening. 
''Observed wavelength of HaA6563. 
'^Observed wavelength of [N II]A6583. 
Emission line values as determined by model fit. 



Note. — The parenthetical velocities are adopted values based on the average FWHM of the other, stronger emission lines. 
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Table 6. CO Emission Line Data 



Source Site line i^v=o VoHsct'' Avks Trms*" ^coAv L'^^" M(H2)'^ 

(GHz) (km s-l) (km s-i) (mK) (Jy km g-^) (xlO^/j-^ (xlO^O/i-^ 

K km pc2) Mo) 



IRAS F09105+4108 


15m 


CO(3 ^ 2) 


239.853 


-104 


37 


1.7 


<7.7 


<3.0 


<1.2 




30m 


CO(3 ^ 2) 


239.853 





31 


6.2 


<8.2 


<3.1 


<1.2 










625 


30 


3.3 


<4.3 


<1.6 


<0.66 




30m 


CO(2 1) 


159.907 





38 


4.6 


<6.4 


<5.6 


<2.2 










937 


38 


3.3 


<4.7 


<4.1 


<1.6 


IRAS F15307+3252 


15m 


CO (4 -» 3) 


239.415 


-108 


38 


1.2 


<5.6 


<4.5 


<1.8 




30m 


CO (4 -» 3) 


239.415 





30 


5.7 


<7.4 


<6.0 


<2.4'' 










188 


30 


8.4 


<11 


<8.9 


<3.5 










-188 


30 


8.4 


<11 


<8.8 


<3.5 










814 


30 


3.6 


<4.7 


<3.8 


<1.5 










-814 


30 


3.0 


<3.9 


<3.1 


<1.2 


PG 1634+706 


12m 


CO(2 1) 


98.647 





49 


0.46 


<3.8 


<23 


<9.2f 


PG 1634+706 


15m 


CO(5 -» 4) 


246.584 





36 


1.6 


<7.4 


<7.1 


<2.8 



^The difference between the velocity the receivers were tuned to and the the systemic velocity of the source. 
''The rms temperatures are given in terms of main beam brightness temperature. 

'^Assuming go = 0.5, Ho = 100 h km s^-*^ Mpc^^. 
■^Assuming a = 4 Mq (K km s"'^ pc'^)"^ (see text). 

"The 3<T upper limit of M(H2) reported by Yun, Scoville, & Evans (1998) over the redshift range 0.9240-0.09275 is 9.8 x lO'^h'^ Mq. 
*The 3(7 upper limit of M(H2) reported by Barvainis et al. (1998) is 2.0 x W^°h-^ Mq. 
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Table 7. Observed Emission- Line Ratios 



Q [OIII]A5007 [QUI] [NII]AA6548+6583 [SII]A6724 



HyLIGs 

IRAS F09105+4108 3.9^ ■•■ 0.55 ± 0.2 0.23 ±0.11 

IRAS F15307+3252 3.0 ± 0.4 > 11.1 ±0.3 1.1 ±0.2 0.19 ±0.10 

IRAS Bright Galaxies 

Seyfert2 1.8 ±1.0 10.2 ± 4.6 1.3 ±0.7 0.41 ± 0.26 

LINERs 0.16 ±0.10 1.7 ±1.0 1.7 ±1.1 0.69 ± 0.41 

HII 0.10 ±0.11 0.64 ±0.52 0.62 ± 0.24 0.27 ± 0.06 



^Determined from [O III] A5007 / H/? = 12 (Kleinmann et al. 1988) and Ha 
/ H/? ~ 3.1. The Ha / H/? ratio is based on the determination by Kleinmann et 
al (1988) that E{B -V)< 0.15 for F09105±4108. 
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Figure Captions 

Figure 1. Broad-band images of the IRAS F09105+4108; a & b) I-band images with 
56" X 56" and ~ 180" x 180" fields of view, respectively, c & d) K'-band images with 
56" X 56" and 143" x 143" fields of view, respectively. For all images. North is up and E is 
to the left. 

Figure 2. Broad-band images of the IRAS F15307+3252; a & b) R-band images with 
60" X 60" and ~ 180" x 180" fields of view, respectively, c & d) K'-band images with 
60" X 60" and 180" x 180" fields of view, respectively. For all images. North is up and E is 
to the left. 

Figure 3. I-band Image of PG 1634+706 with a 72" x 72" field of view. 

Figure 4. KSPEC spectroscopy of IRAS F09105-I-4108 covering the observer-frame 
wavelength range 7000-17500 A. The pixel sampling for the R, 12, Jl, and H-bands are 7.4, 
8.5, 12, and 14 /im pixel"^, respectively. 

Figure 5. KSPEC spectroscopy of IRAS F15307-I-3252 covering the observer-frame 
wavelength range 8600-22400 A. Note that the feature in panel (a) at 17620 A is an OH 
sky line. The pixel sampling for the 12, J2, H, and K'-bands are 8.5, 12, 14, and 20 /xm 
pixer^ respectively. Panel (e) is a CGS4 spectrum of the Ha + [N II] and [S II] A6724 
fines. The pixel sampfing for the CGS4 spectrum is 4.4 /xm pixel"^. 

Figure 6. KSPEC spectrum of [O HI] A5007 and Ha emission from PG 1634+706. The 
pixel sampfing for the H and K'-bands are 12 and 14 //m pixel"^, respectively. Note that 
the broad Fe II complex is also evident redward of 12000A. 

Figure 7. JCMT 15m and IRAM 30m spectra of IRAS F09105+4108 {z - 0.4417). The 
intensity scale is main-beam brightness temperature (IRAM) and antenna temperature 
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corrected for aperture losses (JCMT). A linear baseline has been subtracted from each 
spectrum. The zero velocity corresponds to the average redshift of optical/near- infrared 
emission lines observed for IRAS F09105+4108. 

Figure 8. JCMT 15m and IRAM 30m spectra of IRAS F15307+3252 {z = 0.926). The 
intensity scale is main beam brightness temperature (IRAM) and antenna temperature 
corrected for aperture losses (JCMT). A linear basehne has been subtracted from each 
spectrum. The zero velocity corresponds to the average redshift of optical/near- infrared 
emission lines observed for IRAS F15307+3252. 

Figure 9. NRAO 12m and JCMT 15m spectra of PC 1634+706 {z = 1.338). The intensity 
scale is main beam brightness temperature (NRAO) and antenna temperature corrected for 
aperture losses (JCMT). A linear baseline has been subtracted from each spectrum. The 
zero velocity corresponds to the systemic velocity of the QSO as determined by Tytler & 
Fan (1992) averaged with the redshift of the [O III] A5007 and Ha emission hues. 

Figure 10. Four emission-line diagnostic diagram plots of IRAS F09105+4108 and 
IRAS F15307+3252. The [O I] A6300 / Ha ratio for F09105+4108 has been determined 
from the data in Kleinmann et al. (1988), where the Ha flux was determined from the H/3 
flux and the observation that the extinction in F09105+4108 is neglegible. Additional data 
of H II regions, LINERS, and Seyfert are taken from data compiled by Osterbrock et al. 

(1992) . 

Fi gure 11. (a). L'qq vs. Lir plot of IRAS galaxies and QSOs. Arrows denote upper limits. 
References for CO and infrared (ir) data: LlGs - Sanders, Scoville, & Soifer (1991); ULlGs 
- Sanders (1991); IRAS F09105+4108 - This work (CO) and Kleinmann et al. (1988) (ir); 
IRAS F15307+3252 - This work (CO) and Cutri et al (1994) (ir); IZwl - Barvainis, AUoin, 
& Antonucci (1989); Mrkl014 - Sanders, Scoville, & Soifer (1988a); 3C48 - Scoville et al. 

(1993) ; PC 1634+706 - This work (CO) and Sanders et al. (1989) (ir). (b). LijL'cQ vs. Li, 
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plot of the same sources. Arrows denote upper limits. 

Figure 12. Rest-frame SEDs of IRAS F09105+4108, IRAS F15307+3252, and PG 1634+706. 
The dashed hues represent smooth fits to the data. References for data: For 
IRAS F09105+4108 - all data are from Kleinmann et al. 1988. For IRAS F15307+3252 - all 
data are from Cutri et al. 1994. For PG 1634+706 - all data from Sanders et al. (1989), 
except for the 560 //m data (Chini, Kreysa, & Bierman 1989) and data between 10^-10^//m 
(Antonucci & Barvainis 1988). 
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